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Abstract 

We perform next-to-leading order calculations of the single-diffractive and non-diffractive cross 
sections for dijet production in proton-antiproton collisions at the Tevatron. By comparing their 
ratio to the data published by the CDF collaboration for two different center-of-mass energies, 
we deduce the rapidity-gap survival probability as a function of the momentum fraction of the 
parton in the antiproton. Assuming Regge factorization, this probability can be interpreted as 
a suppression factor for the diffractive structure function measured in deep-inelastic scattering 
at HERA. In contrast to the observations for photoproduction, the suppression factor in proton- 
antiproton collisions depends on the momentum fraction of the parton in the Pomeron even at 
next-to-leading order. 
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I. INTRODUCTION 



Diffractive events in high-energy pp or ep coUisions are characterized by the presence of 
a leading proton or antiproton, which remains intact, and/or by a rapidity gap, defined 
as a (pseudo-)rapidity region devoid of particles. Theoretically, diffractive interactions are 
described in the framework of Regge theory [l| as the exchange of a trajectory with vacuum 
quantum numbers, the so-called Pomeron (iP) trajectory. Diffractive scattering involving 
hard processes (hard diffraction) such as the production of high-i^y jets has been studied 
experimentally to investigate the parton content of the Pomeron (or additional lower-lying 
Regge poles). In this framework, pp hard diffraction can be expressed as a two-step process, 
p + p^p + F + p' —>■ 2 jets + p' + X, and similarly diffractive deep-inelastic scattering 
(DDIS) as Y + P ^ 1* + ^ + P' ^ p' + X. The subprocess 7* iP ^ X is interpreted 
as deep-inelastic scattering (DIS) on the Pomeron target for the case that the virtuality 
of the exchanged photon is sufficiently large. In analogy to DIS on a proton target, 
'~f* + p ^ X, the cross section for DIS on a Pomeron target is expressed as a convolution of 
partonic cross section and universal parton distribution functions (PDFs) of the Pomeron. 
The partonic cross sections are the same as for 7*p DIS. The Pomeron PDFs are multiplied 
with vertex functions for the vertex p ^ IP + p', yielding the diffractive PDFs (DPDFs). 
The additional vertex functions depend on the fractional momentum loss C, and the four- 
momentum transfer squared t of the recoiling proton. The DDIS experiments measure the 
diffractive structure function of the proton F^{^,P,Q^) integrated over t, where /3 = x/^ 
is the momentum fraction of the parton in the Pomeron and is the virtuality of the 7*. 
The evolution of the DPDFs is calculated with the usual DGLAP [2,] evolution equations 
known from 7* — > X DIS. Except for the evolution, the DPDFs are not calculable in 
the framework of perturbative QCD and must be determined from experiment. Such DPDFs 



The presence of a 



have been obtained from the HERA inclusive measurements of yl, 
hard scale such as the squared photon virtuality in DIS or a large transverse jet energy 
E^^^ in hard diffractive processes, as for example in plP jets + X or ^ylP jets + X, allows 
for the calculation of the corresponding partonic cross sections using perturbative QCD. The 
central issue is whether such hard diffractive processes obey QCD factorization, i.e. can be 
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calculated in terms of parton-level cross sections convolved with universal DPDFs. 

For DIS processes, QCD factorization has been proven to hold 15|, and DPDFs have been 

fi n 

extracted at low and intermediate [3|, |J] from high-precision inclusive measurements of 
the process e + p^e'+p' + X using the usual DGLAP evolution equations. The proof 
of the factorization formula also appears to be valid for the direct part of photoproduction 
(Q^ ~ 0) of jets [5j. However, factorization does not hold for hard processes in diffractive 
hadron-hadron scattering. The problem is that soft interactions between the ingoing hadrons 
and their remnants occur in both the initial and final states. This was also the result of 
experimental measurements by the CDF collaboration at the Tevatron [6], where it was 
found that the single-diffractive dijet production cross section was suppressed by up to an 
order of magnitude as compared to the prediction based on DPDFs determined earlier by 
the HI collaboration at HERA p]. In the CDF experiment [gI, the suppression factor was 
determined by comparing single-diffractive (SD) and non-diffractive (ND) events. SD events 
are triggered on a leading antiproton in the Roman pot spectrometer and at least one jet, 
while the ND trigger requires only a jet in the CDF calorimeters. The ratio ^, t) of SD 
to ND dijet production rates A^jj is in a first approximation proportional to the ratio of the 
corresponding structure functions Fjj, i.e. 

An approximation to the SD structure function FjP{x,Q^,^,t), was obtained by 

multiplying the above ratio of rates by the known effective 

F^ri^)=49{x) + ^T.^^i^)] (2) 

after integrating this ratio over ^ and t and changing variables from x to j3 using x — > f3^. The 
result was then compared to the DPDFs from HI 7| using the same approximate formula, 
Eq. (2), relating the structure function to gluon and quark DPDFs as in the ND case. The 
above formula for the ratio R{x, ^, t) is certainly not sufficient for estimating the suppression 
factor for diffractive dijet production in pp collisions. It is based on a leading order (LO) 
calculation of the cross section in the numerator and in the denominator. Furthermore, it 
is assumed that the convolutions of the PDFs in the numerator and the denominator with 
the partonic cross sections are identical and drop out in the ratio together with the PDFs 
for the ingoing proton. These approximations are not valid in next-to- leading order (NLO), 



where, in particular, the cross sections in the numerator and denominator depend on the jet 
algorithm and on the kinematics of the SD and ND processes. 

Since 2002, the two HERA collaborations have presented results for diffractive dijet pho- 
toproduction in order to establish a possible suppression factor. The factorization breaking 
was first investigated on the basis of NLO predictions by us in 2004 j^, by comparing 
to preliminary HI data 10|]. Already in 2004 it became clear that in photoproduction the 
breaking could be shown only by comparing with NLO predictions, which produced by a fac- 
tor of two larger cross sections than the LO predictions. Concerning factorization breaking, 



the conclusions were the same based on a preliminary ZEUS analysis H 
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tions, HI and ZEUS, have now published their final experimental data 
HI confirm in [l^ their earlier findings based on the analysis of their preliminary data and 
preliminary DPDFs, the ZEUS collaboration [l3| reached somewhat different conclusions 
from their analysis. In particular, the HI collaboration 1^ obtained a global suppression of 
their measured cross sections as compared to the NLO calculation of approximately S = 0.5. 
In addition they concluded that also the direct cross section together with the resolved one 
does not obey factorization. The ZEUS collaboration, however, concluded from their analy- 
sis [l^ that, within the large uncertainties of the NLO calculations and the differences in the 
DPDF input, their data are compatible with the NLO QCD calculation, i.e. a suppression 
could not be deduced from their data. In several recent reviews, we have shown, however, 
that the ZEUS data are compatible with the older HI [l^ and with even more recent HI 



data |lJ], if one adjusts the ZEUS large rapidity-gap inclusive DIS diffractive data to the 
analogous HI data, which are the basis of the recent HI DPDFs |4|] and which are used to 
predict the diffractive dijet photoproduction cross sections. In these recent reviews [15|] we 
also investigated whether the NLO prediction with resolved suppression only, which would 
be more in line with the findings in jsj], will also describe the HI and ZEUS data in a satis- 
factory way. The result is, that this is indeed possible, and the resolved suppression factor 
is of the order of S* ~ 0.3. For the global suppression, i.e. direct and resolved component 
equally, the suppression factor is larger, and in addition, depends on Ej^^, which is not the 
case for the resolved suppression only. 

In this work we want to bring the theoretical analysis of diffractive dijet production in 
pp collisions to the same level as has been done for diffractive dijet photoproduction, i.e. to 
calculate the cross sections up to NLO and then compare with the CDF data, to establish 
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the suppression factor in the Tevatron energy range. For this purpose we shall calculate the 
ratio R{x,C,,t). For this we need the NLO cross sections for SD and ND with the cuts as 
in the CDF measurements. The outline of the paper is as follows. In Section 2 we shall 
describe shortly the kinematic restrictions for the CDF analysis based on measurements at 
Run I for y/s = 1800 GeV [6] and on measurements at ^/s = 630 GeV and y/s = 1800 
GeV [16] obtained for comparison at two different center-of-mass energies. In this section 
we shall also specify the various inputs for our calculation. Our results and the comparison 
with the CDF data are presented in Sect. 3. The first 1800 GeV data are compared with 
the calculations in Sect. 3.1. The comparative study of the 630 and the new 1800 GeV cross 
sections are presented in Sect. 3.2. An interpretation of the observed suppression factor is 
given in Sect. 3.3. Sect. 4 contains a summary and our conclusions. 



II. KINEMATIC CUTS AND INPUT FOR THE CALCULATIONS 



The data, which we want to compare our NLO calculations with, are published in Ref. 
6] and Ref. [l^. In the first paper Q|, the CDF collaboration measured non-diffr active 
and single- diffract ive dijet cross sections at a center-of-mass energy of ^/s = 1800 GeV 
using Run IC (1995-1996) data. From an inclusive sample of single-diffraction (SD) events, 
pp — > p'X, triggering on a p detected in a forward Roman pot spectrometer, a diffractive 
dijet subsample with transverse energy ^j?* > 7 GeV was selected. In addition to the two 
leading jets, these events contain other lower- i^-p jets. Similarly, a non-diffractive (ND) dijet 
sample was selected. From the and the rapidity rj of the jets, the fraction Xp of the 
momentum of the antiproton carried by the struck parton was calculated, where Xp is given 
by 

Xp = ^y:Ei,e-^\ (3) 
V -5 



The jets were detected and their energy measured by calorimeters covering the pseudora- 
pidity range |?7| < 4.2. The was defined as the sum of the calorimeter E^s within an 
Tj — (f) cone of radius 0.7. The jet energy correction included a subtraction of an average 
underlying event of Et of 0.54 (1.16) GeV for diffractive (non-diffractive) events. The recoil 
antiproton fractional momentum loss ^ and four-momentum transfer squared t were in the 
range 0.035 < ^ < 0.095 and (— t) < 3 GeV^, respectively, which was in the final sample re- 



stricted to (— t) < 1 GeV^. In the second paper [l6|, the study of diffractive dijet events was 



extended to y/s = 630 GeV. These data were compared to new measurements at ^/s = 1800 
GeV in order to test Regge factorization. This study is similar to the previous diffractive 
dijet study in experimental setup and methodology. For the SD sample, the ^-region is the 
same, 0.035 < ^ < 0.095, but (— t) < 0.2 GeV^. Again in the SD sample events with at least 
two jets with Ej^^ > 7 GeV were selected, where again the Ej^^ was defined as the sum of 
the calorimeter Et^s within a cone of 0.7 in rj — (p space. The jet energy correction included 
a subtraction of an average underlying event of 0.5 (0.9) GeV for SD (ND) events. 

The calculation of the cross sections for dijet production in non-diffractive and single- 
diffractive processes has been performed up to NLO. For the comparison we have calculated 
these cross sections also in LO. For our calculations, we rely on our work on dijet production 
in the reaction 7 + p — jets + X |l7|, in which we have calculated the cross sections for 
inclusive one-jet and two-jet production up to NLO for both the direct and the resolved 
contribution. The version for the resolved contribution can be used immediately for two- 
jet production in pp collisions by substituting for the photon PDF the antiproton PDF 
(for ND) or the Pomeron PDF (for SD). For the (anti-)proton PDF we have chosen the 
version CTEQ6.6M [l^ for the NLO calculation with Nf = 5 active fiavors. The strong 
coupling constant is calculated from the two-loop formula with = 226 MeV. For 
the calculation in LO we have chosen CTEQ6L1 with as determined from the one-loop 
formula and A^^) = 165 MeV. The diffractive PDFs are taken from the recent HI fits to the 
inclusive diffractive DIS data ^. They are only available at NLO and come in two versions, 
'HI 2006 fit A' and 'HI 2006 fit B'. These differ mostly in the gluon density, which is poorly 
constrained by the inclusive diffractive scattering data, since there is no direct coupling of the 
photon to gluons, so that the gluon density is constrained only through the evolution. The 
'HI 2006 fit A' has a much larger gluon for larger momentum fractions (3 at the starting scale 
of Qo = Vs. 5 GeV than 'fit B', which leads to a larger gluon also for larger scales Q. The 
original fit on the data in 41] is performed with Nf = 3 massless fiavors. The production 
of charm quarks was treated in the Fixed-Flavor Number Scheme (FFNS) in NLO with 
non-zero charm-quark mass yielding a diffractive Fg. This E2 is contained in the 'HI fit 
2006 A, B' parameterizations and is then converted into a charm PDF. The HI collaboration 
constructed a third set of DPDFs, which is called the 'HI 2007 fit jets' and which is obtained 
through a simultaneous fit to the diffractive DIS inclusive and dijet cross sections \2^. It 
is performed under the assumption that there is no factorization breaking in the diffractive 
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dijet cross sections. Including the diffractive DIS dijet cross section in the analysis leads 
to additional constraints, mostly on the diffractive gluon distribution. On average, the 'HI 
2007 fit jets' is similar to the 'HI 2006 fit B', except for the gluon distribution at large 
momentum fraction and small factorization scale. The DPDFs of HI contain as a factor the 
vertex function f]p/p{^,t), which describes the coupling of the Pomeron to the proton, i.e. 
the proton-proton-Pomeron vertex. This vertex function is parameterized by the Pomeron 
trajectory ajpit) and an additional exponential dependence on t. This function is used 
also for our calculations, as it has been determined by the HI collaboration when fitting 
their data. The normalization factor of this function is included in the Pomeron PDFs. 
Therefore the HI DPDFs are products of the Pomeron flux factors and the Pomeron PDFs. 
These HI DPDFs include also low-mass proton dissociative processes with invariant mass 
My < 1.6 GeV, which increases the inclusive diffractive DIS cross section as compared to 
cross sections with a pure (anti-)proton final state. We have to keep this in mind, when we 
compare to the CDF data, which use a forward Roman pot spectrometer to trigger on the 
final antiproton and therefore have no antiproton dissociative contributions. 



III. RESULTS 



In this section we present our results and compare them to the experimental data obtained 
with ^/s = 1.8 TeV in 61] and to the more recent data with -^i = 0.63 TeV and ^/s = 1.8 
TeV published in 16(]. In this latter publication, the kinematic constraints differ in some 
points from the constraints used in First we compare to the normalized differential cross 
sections da/dEx and da/dr] for non-diffractive and diffractive dijet production. Second, the 
ratio R{xp) of the number of SD dijet events to the number of ND dijets is compared to the 
CDF data. This function R{xp), obtained by integrating the cross sections in the numerator 
of Eq. (1) over ^ and t, is the main result, and from the theoretical and experimental 
distribution ClS db function of Xp the suppression factor R''''P{xp)/R^^^^^{xp) can be deduced 
and can be studied for the three HI DPDFs, 'HI 2006 fit A', 'HI 2006 fit B' 3 and 'HI 



2007 fit jets' 



20l | for the NLO and the LO (which has been done only for the 'HI 2006 fit 



B') calculations. 
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FIG. 1: Normalized average transverse-energy distributions of the non-diffractive (left) and single- 
diffractive (right) dijet cross section at Run I of the Tevatron. The CDF data (points) are compared 
with our predictions at NLO (full) and LO (dotted) and also with a varied cut on the subleading 
jet Et (dot-dashed). The NLO scale uncertainty is shown as a shaded band (color online). 

A. Comparison with 1800 GeV data 



First we have calculated the distribution ^-^= as a function of Et = {E^^^^ + i?:|?*^)/2, 
where E^r^^^ [E^r^^^'^) refers to the jet with the largest (second largest) Et for ND and SD dijet 
production with ^/s = 1800 GeV center-of-mass energy, integrated over the rapidities of 
the jets in the range \7]\ < 4.2. Jets are defined with the usual cone algorithm within a 



chosen r] — cj) cone of radius R = 0.7 and a partonic distance R^ep = 1.3R to match the 
experimental analysis 2l|. a is the integrated cross section with the cut E^r^^^^'^'^ > 7.0(6.5) 
GeV. The lower limit of the leading and subleading jet differ slightly in order to avoid 
infrared sensitivity in the computation of the NLO cross sections, when integrated over Et 



22l |. Unfortunately in the experimental analysis such an asymmetric choice of E^^^^ and 



E}^^'^ has not been made, since both E^^^^ and E^r^^"^ are restricted by E^^^^''^ > 7.0 GeV, so 
that we do not know whether the choice E^t^'^ > 6.5 GeV is in accord with the experimental 
analysis. Therefore we have also varied the E^r^^^ cut slightly to E'^*^ > 6.6 GeV. The results 
for ^-^= are shown in Fig. 1 (left) for E';!?*^ > 6.5 GeV (full histogram) and i?^*^ > 6.6 
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GeV (dot-dashed histogram), respectively. Together with the NLO cross section, we also 
show the LO cross section (dotted histogram) and the scale variation of the NLO result 
(shaded band), where the renormalization and factorization scales are varied simultaneously 
by factors of 0.5 and 2.0 compared to the default scale, which is chosen equal to E^r^^^, i.e. the 
largest Et of both jets. As is often the case, the scale uncertainty is relatively small in the 
normalized distributions. In Fig. 1 (left) we have included also the measured cross section 
from [q, which unfortunately is given without the experimental uncertainties. Besides the 
statistical errors, which should be similar to those in the single-diffractive distributions given 
the similar number of ND and SD events 211], there are also systematic errors, as for example 
those associated with the jet energy scale. These would be needed for a fair comparison. 
Second, the theoretical cross sections should be corrected for hadronization effects when 
comparing to data. These are not known to us, but could be calculated through Monte 
Carlo models. If we compare the calculations in Fig. 1 (left) with the data, we observe 
that the results with > 6.5 GeV agree reasonably well with the data in the large Et 
range, Et > 20 GeV, but much less for the low and medium Et range. Conversely, for 
£"2?*^ > 6.6 GeV the small and medium Et range agrees better and the large Et range less. 
The experimental errors will be larger in the large Et range. Therefore the cross section 
with the cut > 6.6 GeV would be preferred, in particular also because we have perfect 
agreement in the first, second and third Et bin, which are the most important ones for the 
integrated cross section cr, which determines the normalization. 

The equivalent comparison for the SD dijet ii^r-distribution is shown in Fig. 1 (right) for 
E^T^^^^ > 7.0(6.5) GeV (full) and E^^"^^^ > 7.0(6.6) GeV (dot-dashed). Here we have chosen 
only the 'HI 2006 fit B' as DPDF. The comparison of data versus theory in Fig. 1 (right) 
shows the same pattern as for the ND distributions in Fig. 1 (left). In general the agreement 
with the data is even somewhat better now, in particular for the E^^"^^^ > 7.0(6.6) GeV 
cut. As for the ND distribution, we present also the LO prediction (dotted). From the 
unnormalized distributions (not shown), we obtain ratios of NLO to the LO cross sections 
(i^-factors), which increase from relatively small values of 0.5 (0.6) in the infrared-sensitive 
region close to the E^r^^^''^ cuts to unity at larger Et for the ND (SD) cross sections, indicating 
good perturbative stability and no sensitivity to the cut on E^r^^'^ there. 

The ND fy-distribution where fj = (ryj^*^ + ri^'^^^)/2, is plotted in Fig. 2 (left), again 
for the two cuts, and the SD ^-distribution in Fig. 2 (right). For ND and SD, the two 
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FIG. 2: Same as Fig. 1, but for the normalized average rapidity distributions (color online). 

choices for the E^r^^'^ have little influence on the result. The cross sections are somewhat 
smaller (by about 20%) at the maximum as compared to the ND experimental data, again 
given without experimental errors, but are in agreement with the SD data. The theoretical 
diffractive 77-distribution is boosted towards positive f/, as is the experimental one. We 
conclude from these comparisons that there is satisfactory agreement between the measured 
Et and rj distributions and our theoretical predictions based on the 'HI 2006 fit B' DPDF. 
This motivates us to move on to the calculation of the ratio R{xp) of SD to ND cross sections. 
Actually, if the experimental cross sections above had been known to us in the unnormalized 
form, we would have been in the position to deduce suppression factors as a function of inl- 
and f]. 

The ratio R{xp) of the SD to ND cross sections is evaluated as a function of Xp, the fraction 
of the momentum of the antiproton transferred to the struck parton. It is calculated from 
the E^ and r/* of the jets with the relation in Eq. (3), where the sum over i is taken over 
the two leading jets plus the next highest Et jet with Et > 5 GeV. The cross section 
in the numerator is integrated over ^ in the range 0.035 < < 0.095 and over (— t) in 
the range (— t) < 1 GeV^ and over the Ej^^ of the highest and second highest Ej^^ with 
^jeti(2) ^ 7.0(6.5) GeV for both the SD and ND jet sample. As already mentioned, the 
data samples have the constraint E^r^^^^"^^ > 7.0(7.0) GeV. We have checked that the choice 
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FIG. 3: Left: Ratio K of SD to ND dijet cross sections as a function of the momentum fraction of the 
parton in the antiproton, computed at NLO (with three different DPDFs) and at LO and compared 
to the Tevatron Run I data from the CDF collaboration. Right: Double ratio of experimental over 
theoretical values of i?, equivalent to the factorization-breaking suppression factor required for an 
accurate theoretical description of the data (color online). 

Ef^ > 6.6 GeV for the ND and SD cross sections has negligible influence on R{xp). The 

results are plotted in Fig. 3 (left) as a function of log^gl^p) three choices of the DPDFs, 

'HI 2006 fit B' (full), 'HI 2006 fit A' (dashed) and 'HI 2007 fit jets' (dot-dashed). All three 

are calculated in NLO. The NLO scale uncertainty for 'HI 2006 fit B' (shaded band) cancels 

out to a large extent in this ratio of cross sections. The LO prediction for 'HI 2006 fit B' 

fl 

is also given (dotted). The CDF data, which are plotted in Ref. [g] in six ^ bins of width 
A,^ = 0.01, have been integrated to give R{xp) in the range 0.035 < ^ < 0.095. They were 
available in numerical form with statistical errors 23|] and are also plotted in Fig. 3 (left). 
From these presentations it is obvious that the theoretical ratios R{xp) are, depending on 
Xp, by up to an order of magnitude larger than the measured R{xp) in agreement with the 
result in 61]. There are quite some differences for the different DPDF choices. In general 
'fit B' and 'fit jets' lie closely together, whereas 'fit A' deviates more or less from these two 
depending on the Xp range. For 'fit B' we also show the scale error and the LO prediction. 
The hierarchy between the three DPDFs at large Xp is easily explained by the fact that the 



11 



corresponding gluon DPDFs are at large Xp the largest for 'fit A' and the smallest for 'fit 
jets' 0, . The same pattern between the different DPDFs is seen even more clearly if we 
plot the ratio of the experimental R{xp) and the theoretical R{xp) as a function of log^o(^p)- 
The result for this (double) ratio R'^^^{xp)/ R^^^^'^{xp) is seen in Fig. 3 (right). As can be 
seen, this ratio varies in a rather similar way for the three DPDFs in NLO and for 'fit B' in 
LO in the range 10^^ < Xp < 10^^. The variation is strongest for the 'fit A' DPDF, where 
this ratio varies by more than a factor of seven. For the other two DPDFs this variation 
is somewhat less, but still appreciable. Actually, we would expect that the ratio plotted in 
Fig. 3 (right), which gives us the suppression factor, should vary only moderately with Xp. 
After presenting the y/s = 630 GeV and the more recent y/s = 1800 GeV data below, we 
shall discuss possible interpretations of this result. We also observe that the suppression 
factor for 'fit B' in NLO and LO are different, in particular for the very small Xp. 

In Ref. (gI, the ratio R{xp) was multiplied with an effective PDF governing the ND 
cross section to obtain the effective DPDF -Fjj(/?) as a function of P = Xp/^. This effective 
non- diffract ive PDF Fjj^(x) is calculated from the formula in Eq. (2), where the gluon PDF 



g{x) and the quark PDFs qi{x) are taken from the GRV98 LO parton density set 2^ and 
evaluated at the scale = 75 GeV^, corresponding to the average E}^' of the SD and ND 
jet cross sections. Then, the effective diffractive PDF F^j{f3) of the antiproton is obtained 
from the equation 

F°(/?) = R{x = (30 X Fjr(x PO- (4) 

We use this relation for the experimental and theoretical values of R{xp). However, both are 
integrated over ^ and are not given as function of ^. We consider them as only moderately 
^-dependent and evaluate R{x = (30 -Fjj°(a; — > /3,^) at an average value of ^ = 0.0631. 
This works quite well over the /?-range of interest, if we compare the Fjj(/?) values obtained 
in this way with the FPj{(3) in the CDF publication |6|], as is seen in Fig. 4 (left), where 
the full published points from 6j, denoted CDF PRL 84, coincide rather well with the open 
points deduced from the above equation and the published values of R, denoted CDF PRL 84 
(KK). The ratio F°'""P(/5)/F°'^^^^°(/5), plotted in Fig. 4 (right) clS 8b function of (3 linearly, 
gives the suppression factor as a function of (3 instead of Xp cLS 111 Fig. 3. For example, for 
the 'fit B' DPDF it varies between 0.3 at (3 = 0.05 and 0.13 at /3 = 0.1 to 0.07 at (3 = 0.9. In 
the range 0.3 < (3 < 0.9 the suppression factor varies only moderately with /5, but increases 
strongly for (3 < 0.3, independently of the chosen DPDF. Above (3 = 0.3, 'fit B' and 'fit 
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FIG. 4: Left: Effective diffractive structure function i<jj of the partons with momentum frac- 
tion l3 in the Pomeron as measured in dijet production at the Tevatron and compared to our 
(N)LO calculations. Right: Ratio of experimental over theoretical values of -Pjj, equivalent to 
the factorization-breaking suppression factor required for an accurate theoretical description of the 
data (color online). 

jets' show the most constant behavior. Here one should note that the result in Fig. 4 (right) 



is independent of the assumptions inherent in Eq. (4), since 



ND, 
JJ ' 



X 



P^) cancels in the 



ratio. The information in this figure concerning the suppression factor is equivalent to Fig. 
4 of the CDF publication 6|. The main difference to the CDF plot is the fact that now the 
suppression factor is given by comparing to calculated NLO cross sections without using the 
approximate formula Eq. (4) above, which can be justified only in LO. 

To obtain an idea how large the effect of our NLO dijet evaluation compared to a sim- 
ple combination of LO parton densities in the Pomeron is, we have calculated the ratio 
Fj^''^^*~'(/5)/Fj^''"^(/3) for the three DPDFs. Here the numerator is the Fj] from Eq. (4) with 
R evaluated in NLO, i.e. Fig. 4 (left), and the denominator is -Pjj''"^(/?) calculated from the 
formula 



f^D.LO^ 



4 



'{(3) = Jd^J dtfjp/p{^, t)(3[g{(3) + - ^ (5) 
where the Pomeron flux factor fip/p{C,,t) and the gluon and quark PDFs in the Pomeron 
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g{f3) and qi{f3) are taken from the fits 'HI fit A', 'HI fit B' and 'HI fit jets' at the scale 
= 75 GeV^, respectively. At (3 = 0.1, we obtain ratios of 0.95, 1.05 and 1.1 for these three 
fits, respectively, indicating that our more accurate NLO calculations lead to very similar 
suppression factors as the simple approximation in Eq. (5) for all three DPDFs. This ratio 
is more or less constant as a function of (3 in the considered range, meaning that already 
in the CDF publication [6] one has the strong variation of the suppression factor with (3 
as mentioned above. It is interesting to note that replacing the approximate Eq. (5) with 
the experimentally used Eq. (4) compensates the effect of the NLO corrections, as the ratio 
of SD to ND i^"-factors, or equivalently the ratio of the NLO over the LO value of R, is 
approximately 1.35 for the 1800 GeV calculation discussed here and 1.6 for the 630 GeV 
calculation presented in the next subsection. To compute the effect of this approximation 
alone, i.e. the ratio of Eq. (4) at LO over Eq. (5), one must divide the values of 0.95, 1.05 
and 1.1 by the ratio of if-factors, i.e. 1.35. 



B. Comparison with 630 GeV and new 1800 GeV data 

In a second publication, the CDF collaboration presented data for diffractive and non- 
diffractive jet production at ^/s = 630 GeV and compared them with a new measurement at 
^/s = 1800 GeV jl6|. From both measurements they deduced diffractive structure functions 
using the formula Eq. (4) with the expectation that -Fjj(/5) is larger at ^/s = 630 GeV than 
at ^/s = 1800 GeV. The experimental cuts are similar to the cuts in the first analysis 6| 
with the exception that now (— t) < 0.2 GeV^ and in addition to the E^^^^''^ > 7.0 GeV cut 
they require > 10 GeV. This second cut on Erp is very important for the comparison 
with the NLO predictions, since with this additional constraint the infrared sensitivity is 
not present anymore. 

With these cuts and the integration over ^ in the range 0.035 < ^ < 0.095, we have 
calculated the normalized cross sections {l/a)da/dEjn and (l/(T)dcr/dr/ as in the previous 
subsection, but now for y/s = 630 GeV. For ND (left) and SD (right) jet production, the 



results are shown in Fig. [5] and compared to the data from Ref. 16|. Here, the ND data 



sample was larger by about two orders of magnitude compared to the SD data sample, so 



that the statistical errors, which were not given in Ref. 16|, should be smaller by about a 



factor of ten 



2l| . Again, no information about systematic errors was available. We find 
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FIG. 5: Same as Fig. 1, but for a reduced center-of-mass energy of 630 GeV at the Tevatron (color 
online) . 

reasonably good agreement in the medium-i?T range. In these figures, we have also plotted 
the LO predictions (dotted). For the DPDF, we have chosen as before the 'HI 2006 fit B' 
set. Due to the large experimental errors for Ex > 15 GeV for the SD case, we also find 
good agreement in the laige-Ex range. 

The equivalent result and comparison with the data for the rj-distribution is shown in 
Fig. El again for ND (left) and SD (right) jet production for NLO (full) and LO (dotted) 
predictions. The agreement between the theoretical results and the CDF data is similar as 
in the previous subsection, where we compared to the y/s = 1800 GeV data. This justifies 
to go on with the calculation of the ratio R{xp) of SD to ND cross sections. 

The momentum fraction Xp is calculated as before from Eq. (3), and then the cross 
sections da/dxp can be calculated with the same restrictions on the number of included 
jets as before. The only difference is the different cut on (— t). The results for R{xp) at 
y/s = 630 GeV are presented in Fig. [7] (left) for the three choices of DPDFs, 'fit B' (NLO 
and LO), 'fit A' (NLO) and 'fit jets' (NLO). In this figure, also the experimental data from 
Ref. y are included. The range of Xp is now much smaller than for the ^/s = 1800 GeV 
case. It ranges from Xp = 0.025 to Xp = 0.1. From this plot, the suppression of the SD cross 
section is clearly visible. The suppression factor is of the same order of magnitude as in the 
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FIG. 6: Same as Fig. 2, but for a reduced center-of-mass energy of 630 GeV at the Tevatron (color 
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FIG. 8: Double ratios of experimental over theoretical values of R, equivalent to the factorization- 
breaking suppression factor required for an accurate theoretical description of the data from the 
Tevatron at ^/s = 630 (left) and 1800 GeV (right) (color online). 



previous subsection. The same plot for the new i/i = 1800 GeV data [l6| together with the 
predictions is given in Fig. [7] (right). 

From the two plots in Fig. [7] we have calculated the corresponding suppression factors 
^^^P(xp)/^(N)LO(xp), exhibited in Fig. E (left: ^ = 630 GeV; right: ^ = 1800 GeV). 
In both figures we observe that the LO and NLO results for the suppression factors differ 
significantly (LO only given for 'fit B'), but also the three different DPDFs give different 
suppression factors, although with smaller variation compared to the LO and NLO result. 
Due to the variation of this factor with Xp it is difficult to compare the suppression of the 
^/s = 630 GeV result (left) with the i/i = 1800 GeV result (right) in Fig. [HI On average, 
it seems that for larger Xp the two suppression factors are more or less equal and we cannot 
say that the suppression factor for y/s = 630 GeV is larger than for y/s = 1800 GeV, as 
we would expect it. In the region Xp > 0.02, the suppression factors for both y/s are fairly 
constant (~ 0.05), in particular for the DPDF 'fit jets'. This is not the case for the analysis 
in the previous subsection, where, as we see in Fig. 4, the suppression factor varies already 
much more in this particular 

From the results in Fig. [71 we have calculated Fj]{j3) by changing variables from Xp to (3 
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FIG. 9: Effective difFractive structure function Fjj of the partons with momentum fraction /? in 
the Pomeron as measured in dijet production at the Tevatron with y/s = 630 (left) and 1800 GeV 
(right) and compared to our (N)LO calculations (color online). 



with (3 = Xp/^ and = 0.0631 and multiplying with the effective PDF for ND jet production 
as in Eq. (4). The results, together with the corresponding experimental data from Ref. 16 1 
and those calculated with the chosen ^, which agree inside errors except for two points at 
small are shown in Fig. M (left: y/s = 630 GeV, right: y/s = 1800 GeV). 

From these results we have again calculated, as in the previous subsection, the suppression 
factor as a function of (3 in the range < (3 < 0.8. The plots for the ratios Fjj ^/f}^^'"'^ are 
seen in Fig. [10] for the lower (left) and the higher center-of-mass energy data (right), again 
for the three DPDF fits in NLO and 'fit B' also in LO. First we observe that the ratios in 
Fig. [TD] differ very little, except perhaps at very small (3. This means that from these data 
there is no essential difference seen in the suppression at y/s = 630 GeV and y/s = 1800 
GeV. Second, we notice that with the 'fit jets' we have the most constant behavior of the 
suppression for f3 > 0.2. Furthermore, comparing Fig. [10] (right) with Fig. 4 (right) we see 
some differences. While the general pattern is the same, the suppression factor for 'fit jets' 
in particular is less constant and larger in Fig. 4 (right) than in Fig. [10] (right), which is 
obviously correlated with the more restrictive cuts on Ej- and t in the latter. 

For completeness we also compared our NLO dijet calculation to the approximate LO 
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FIG. 10: Ratios of experimental over theoretical values of Fjj for ^/s = 630 (left) and 1800 
GeV (right), equivalent to the factorization-breaking suppression factors required for an accurate 
theoretical description of the data (color online) . 

formula in Eq. (5). For y/s = 1800 GeV and (3 = 0.1, we obtain the same values of 0.95, 
1.05, and 1.1 for 'HI 2006 fit A', 'HI 2006 fit B', and 'HI 2007 fit jets' as for the older CDF 
analysis. They depend again weakly on f3. For ^/s = 630 GeV and /? = 0.1 we obtain larger 
values of 1.15, 1.35, and 1.45, which is in line with the larger ratio of fT-factors (1.6 instead 
of 1.35) for SD and ND events at this lower center-of-mass energy. 

As stated above, the calculation of the effective diffractive structure function Fjj{f3) from 
the ratio R{xp) was based on the assumption that the latter was only weakly ^-dependent, so 
that Eq. (4) could be evaluated at an average value of ^ = 0.0631. This weak ^-dependence 
is indeed observed in the newer CDF data, published in the lower part of Fig. 4 of Ref. [l^ 
and reproduced in our Fig. [H] (full circles). These data agree well with the ^-dependent 
values of R{xp) published in Fig. 3 of Ref. [6[ when transformed into Fjj{^) using Eq. (4) and 
C = Xp/P with P = 0.1 (open circles). The same weak ^-dependence is also observed in our 
theoretical calculations when using the same procedure, except with different normalization, 
reflecting the ^-dependence of the HI fits to the Pomeron flux factors fip/p{C,,t) oc with 
m ~ 1.1 (0.9 in the CDF fit to their data). At the considered value of /? = 0.1, the NLO 
suppression factors for 'fit A,B' and 'fit jets' are 0.15, 0.12 and 0.11, respectively, and are 
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FIG. 11: Effective diffractive structure function in Eq. (4) from (N)LO dijet cross sections for fixed 
/3 as a function of the momentum fraction of the Pomeron in the antiproton ^, compared to the 
Tevatron data at -^s = 1800 GeV from the CDF collaboration 161] (color online). 



almost independent of ^. At LO, the suppression factor for 'fit B' is larger, i.e. 0.15, which 
corresponds to the fact that the ratio of SD over ND i^-factors is 1.35. Note that Fig. [TTlis 
based on the higher statistics CDF data without the stronger cuts on Et and t and should 
therefore be not compared to Fig. [10] (right), but to Fig. 4 (right), where consistency of 
these numbers with the values shown at /5 = 0.1 can be found. 
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The (small) difference of the theoretical (1-1) and experimental (0.9) values of m can 
be explained by a subleading Reggeon contribution, which has not been included in our 
predictions. To study its importance, we have computed the ratio of the Reggeon over 
the Pomeron contribution to the LO single-diffractive cross section at ^/s = 1800 GeV as 
a function of Et, f], and Xp. The Reggeon flux factor was obtained from 'HI 2006 fit B' 
and convolved, as it was done in this fit, with the parton densities in the pion of Owens 



251]. Very similar results were obtained for the 'HI 2007 fit jets' Reggeon flux. On average, 



the Reggeon adds a 5% contribution to the single-diffractive cross section, which is almost 
independent of Et, but is smaller in the proton direction at large f] and small Xp = ^j3 and 
^ (2.5%) than at large Xp and ^ (8%). This corresponds to the graphs shown in Figs. 5 
= 0.01) and 6 (^ = 0.03) of the HI pubhcation [j, e.g. at = 90 GeV^. While the 
Reggeon contribution thus increases the diffractive cross section and reduces the suppression 
factor at large Xp in Figs. 3, 8, and 9, making the latter more constant, the same is less true 
at small values of a;^. 



C. Interpretation of the observed suppression factor 



Our main results are the plots for the suppression factors as a function o f log mfxp) in Fig. 
3 (right) deduced from the data of Ref. [6i] and in Fig. [Slfrom the data of [16'| at ^/s = 630 
GeV (left) and ^/s = 1800 GeV (right). The qualitative behavior of the suppression factor in 
these three figures is very similar. We observe an appreciable dependence of the suppression 
factor on the chosen DPDFs and a dependence on Xp with a minimum at Xp ~ 0.032 
(log]^Q(xp) ~ —1.5) and a rise towards smaller Xp by up to a factor of five. The equivalent 
result as a function of (3 is shown in Fig. 4 (right) for the data of Ref. ^ and in Fig. [10] 
(left) for the = 630 GeV data of Ref. (lej and in Fig. M (right) for the = 1800 GeV 
data of Ref. 16[. Depending on the chosen DPDFs, the suppression factor as a function 
of f3 is minimal with the value ~ 0.05 at f3 = 0.5 and rises with decreasing /5 to a value 
~ 0.12 in Fig. 4 (right) and Fig. [10] (right) and to a value ~ 0.1 in Fig. [10] (left) at (3 = 0.1 
(considering 'fit B' as an example). Of course, this rise of the suppression factor towards 
small (3 is directly related to its rise as a function of Xp towards small Xp. 

A comparison of the HI data 4], which are used to obtain the DPDFs applied in our 
calculation, with a similar measurement, in which the leading proton is directly detected 
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261], yields a ratio of cross sections for My < 1.6 GeV and My = rup of 1.23 ± 0.03{stat.) ± 
0.1Q{syst.) {4]. Since the CDF measurements are performed by triggering on the leading 
antiproton, these measurements must be multiplied by this ratio to normalize them to the 
My < 1.6 GeV constraint for the HI DPDFs. Therefore, all suppression factors obtained so 
far must be multiplied by this ratio. 

Any model calculation of the suppression factor, which is also sometimes called the rapid- 
ity gap survival factor, must try to explain two points, first the amount of suppression, which 
is ~ 0.1 at /9 = 0.1, and second its dependence on the variable P (or Xp). Such a calculation 
has been performed by Kaidalov et al. 27|]. In this calculation, which we call KKMR, the 
hard scattering cross section for the diffractive production of dijets was supplemented by 
screening or absorptive corrections on the basis of eikonal corrections in impact parameter 
(6) space. The parameters of the eikonal were obtained from a two-channel description of 
high-energy inelastic diffraction. The exponentiation of the eikonal stands for the exchange 
of multi-Pomeron contributions, which violate Regge and QCD factorization and modify 
the predictions based on single Pomeron and/or Regge exchange. The obtained suppression 
factor S is not universal, but depends on the details of the hard subprocess as well as on 
the kinematic configurations. The first important observation in the KKMR analysis is that 
in the Tevatron dijet analysis the mass squared of the produced dijet system Mjj = XpP^s 
as well as ^ are almost constant, so that small (3 implies large Xp. The second important 
ingredient in the KKMR model is the assumption that the absorption cross section of the 
valence and the sea components, where the latter includes the gluon, of the incoming pro- 
ton are different, in particular, that the valence and sea components correspond to smaller 
and larger absorption. For large Xp or small f3, the valence quark contribution dominates, 
which produces smaller absorptive cross sections as compared to the sea quark and gluon 
contributions, which dominate at small Xp. Hence the survival probability (or suppression 
factor) increases as Xp increases and P decreases. In Ref. [27|, the convolution of the old 
HI DPDFs [281] and the /3-dependent absorption corrections produced a Fj^(/5)-distribution 
corrected for the soft rescattering, which was in very good agreement with the corresponding 
experimental distribution in the CDF publication |6| (see Fig. 4 in 27]). We have no doubt 



that using our single-diffractive NLO cross sections based on the more recent DPDFs of HI 
[4] will lead to a very similar result. 

An alternative model for the calculation of the suppression factor was developed by 
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Gotsman et al. [29|. However, these authors did not convolve their suppression mechanism 
with the hard scattering cross section. Therefore a direct comparison to the CDF data is 
not possible. 

At variance with the above discussion of diffractive dijet production in hadron-hadron 
scattering, the survival probability in diffractive dijet photoproduction was found to be 
larger (~ 0.5 for global suppression, ~ 0.3 for resolved photon suppression only) and fairly 



independent of (3 (or zjp) 12|, ll5j. This can be explained by the fact that the HERA analyses 



are restricted to large values of > 0.1 (as opposed to small and intermediate values of Xp = 
0.02 ... 0.2), where direct photons or their fluctuations into perturbative or vector meson- 
like valence quarks dominate. The larger suppression factor in photoproduction corresponds 
also to the smaller center-of-mass energy available at HERA. 



IV. CONCLUSIONS 



In conclusion, we have performed the first next-to-leading order calculation of single- 
diffractive and non-diffractive cross sections for dijet production in proton-antiproton col- 
lisions at the Tevatron, using recently obtained parton densities in the (anti-)proton from 
global fits and in the Pomeron from inclusive deep-inelastic scattering and DIS dijet produc- 
tion at HERA. The normalized distributions in the average transverse energy and rapidity 
of the two jets agreed well with those measured by the CDF collaboration at two different 
center-of-mass energies of ^/s = 1800 and 630 GeV. However, the ratios of single-diffractive 
and non-diffractive cross sections had two be multiplied by factors of about 0.05 and up to 
0.3, depending on the momentum fraction of the parton in the antiproton, the center-of- 
mass energy, the order of the calculation, and the DPDF. Assuming Regge factorization, the 
ratios of cross sections were interpreted as ratios of effective diffractive structure functions, 
exhibiting similar suppression factors. 

We found that the ratios of SD over ND i^'-factors of 1.35 and 1.6 at ^/s = 1800 and 630 
GeV, respectively, were partially compensated by the simplification inherent in the definition 
of the effective structure functions, but that the suppression factors were still smaller at NLO 
than at LO. They were also less dependent on the momentum fraction of the parton in the 
Pomeron at NLO than at LO, in particular at the lower center-of-mass energy and to a 
smaller extent also for the more restricted kinematics at the higher ^/s. The DPDF fit 
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by the HI collaboration using DIS dijet data to better constrain the gluon density in the 
Pomeron showed the most constant behavior. 

We pointed out that all suppression factors obtained so far must be corrected by a factor 
of 1.23 ± 0.03{stat.) ± 0.16{syst.) due to the fact that the DPDFs were obtained from HI 
data that includes diffractive dissociation, while the CDF data were triggered on a leading 
antiproton. We also recalled that the remaining momentum-fraction dependence can be 
explained by a two-channel eikonal model that predicts different behaviors for the regions 
dominated by valence quarks and sea quarks and gluons in the proton. This is in contrast 
to the constant behavior observed in photoproduction, which is governed by direct photon 
or valence-like quark contributions. Wc finally confirmed that the single-diffractive data are 
dominated by a single Pomeron exchange, since its momentum fraction dependence in the 
antiproton is well described in shape by the Pomeron flux factors fltted to the HI DIS data. 
An additional Reggeon exchange would add only two to eight percent to the single-diffractive 
cross section for smaller and larger values of the Reggeon momentum fraction. 

Acknowledgment s 

We thank K. Hatakeyama for useful discussions concerning the CDF data analyses. This 
work has been supported by the Theory-LHC- France initiative of the CNRS/IN2P3. 



[1] P.D.B. Collins, An Introduction to Regge Theory and High-Energy Physics, Cambridge Uni- 
versity Press, Cambridge (1977). 

[2] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 5, 438 (1972) and Sov. J. Nucl. Phys. 20, 
94 (1975); G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977); Y.L. Dokshitzer, Sov. 
Phys. JETP 46, 641 (1977). 

[3] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 38, 43 (2004). 

[4] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 48, 715 (2006). 

[5] J.C. Collins, Phys. Rev. D 57, 3051 (1998) [Erratum ibid. 61, 019902 (2000)]; J. Phys. G 28, 
1069 (2002). 

[6] CDF Collaboration, T. Affolder et al., Phys. Rev. Lett. 84, 5043 (2000). 



24 



[7] HI Collaboration, paper 980, submitted to the 31st Int. Conf. on High-Energy Physics 

(ICHEP02), Amsterdam (2002). 
[8] M. Klasen and G. Kramer, in Proc. of the 12th Int. Workshop on Deep Inelastic Scattering 
(DIS04), eds. D. Bruncko, I. Ferencei and P. Strizenec, Kosize, Inst. Exp. Phys., SAS, p. 492 
(2004), |hep-ph/0401202 
[9] M. Klasen and G. Kramer, Eur. Phys. J. C 38, 93 (2004). 

[10] HI Collaboration, paper 987, submitted to the 31st Int. Conf. on High-Energy Physics 
(ICHEP02), Amsterdam (2002); paper 087, submitted to the Int. Europhysics Conf. on High- 
Energy Physics (EPS03), Aachen (2003); F.P. Schilling, Eur. Phys. J. C 33, S530 (2004). 

[11] ZEUS Collaboration, abstract 6-0249, contributed to the 32nd Int. Conf. on High-Energy 
Physics (ICHEP04), Beijing (2004). 

[12] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 51, 549 (2007). 

[13] ZEUS Collaboration, S. Chekanov et al., Eur. Phys. J. C 55, 177 (2008). 

[14] HI Collaboration, K. Cerny et al. in Proc. of the 16th Int. Workshop on Deep-Inelastic 
Scattering (DIS08), London (2008), [http : //dx . doi . org/10 . 3360/dis . 2008769 

[15] M. Klasen and G. Kramer, Mod. Phys. Lett. A 23, 1885 (2008); Proc. of the Workshop on 
HERA and the LHC, eds. H. Jung et al., DESY-PROC-2009-02, Mar 2009. larX iv:0903.386H 
larXiv:0808.3700l LPSC 08-115 (to be published); G. Kramer, Nucl. Phys. B (Proc. Suppl.) 
191, 231 (2009). 

[16] CDF Collaboration, D. Acosta et al., Phys. Rev. Lett. 88, 151802 (2002). 

[17] M. Klasen, T. Kleinwort and G. Kramer, Eur. Phys. J. direct C 1, 1(1998) and the earlier 

papers quoted there; for a review see M. Klasen, Rev. Mod. Phys. 74, 1221 (2002). 
[18] CTEQ Cohaboration, P.M. Nadolsky et al., Phys. Rev. D 78, 013004 (2008). 
[19] CTEQ Collaboration, I. Pumplin et al., JHEP 07, 021 (2002). 
[20] HI Collaboration, A. Aktas et al., JHEP 10, 042 (2007). 
[21] K. Hatakeyama, private communication. 

[22] M. Klasen and G. Kramer, Phys. Lett. B 366, 385 (1996); S. Frixione and G. Ridolfi, Nucl. 

Phys. B 507, 315 (1997). 
[23] See http : /physics . rockefeller . edu/hatake/phys/ sdj j_1800_prl . html. 
[24] M. Gliick, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 (1998). 
[25] J. Owens, Phys. Rev. D 30, 943 (1984). 



25 



[26] HI Collaboration, A. Aktas et al., Eur. Phys. J. C 48, 749 (2006). 

[27] A.B. Kaidalov, V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C 21, 521 (2001). 
[28] HI Collaboration, T. Ahmed et al., Phys. Lett. B 348, 681 (1995); HI Collaboration, C. 

AdlofFet al., Z. Phys. C 76, 613 (1997). 
[29] G. Gotsman, E. Levin, U. Maor, E. Naftali and A. Prygarin, Proc. of the Workshop on HERA 

and the LHC, part A, p. 221 (2005). 



26 



